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Abstract 

In this paper we derive the canonical distribution as a stationary solution of the Liouville equation 
for the classical dissipative system. Dissipative classical systems can have stationary states look like 
canonical Gibbs distributions. The condition for non-potential forces which leads to this stationary 
solution is very simple: the power of the non-potential forces must be directly proportional to the 
velocity of the Gibbs phase (phase entropy density) change. The example of the canonical distribution 
for a linear oscillator with friction is considered. 
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1 Introduction 

The canonical distribution was defined 101 years ago in the book "Elementary principles in statistical 
mechanics, developed with especial reference to the rational foundation of thermodynamics" [1] published 
in 1902. The canonical distribution function usually can be derived as a stationary solution of the Liouville 
equation for non-dissipative Hamiltonian N-particle systems with a special set of potential forces [2]- [6]. 

In general, classical systems are not Hamiltonian systems and the forces which act on particles are 
the sum of potential and non-potential forces. The non-potential internal forces for an N-particle system 
can be connected with nonelastic collisions [7]. Dissipative and non-Hamiltonian systems can have the 
same stationary states as Hamiltonian systems [8]. For example, dissipative quantum systems have pure 
stationary states of linear harmonic oscillator [9]. We can assume that canonical distribution exists for 
classical dissipative systems. 

In this paper we consider the Liouville equation for classical dissipative and non-Hamiltonian N- 
particle systems. This equation is the equation of continuity in 6N dimensional phase space. We find 
the condition for the non-potential forces which leads to the stationary solution of this equation look like 
the canonical distribution function. This condition is very simple: the power of non-potential forces must 
be directly proportional to the velocity of the Gibbs phase (phase entropy density) change. Note that the 
velocity of the phase entropy density change is equal to the velocity of the phase volume change. 

In Sections 2, the mathematical background and notations are considered. In this section, we formu- 
late the main conditions for non-potential forces and derive the canonical distribution from an N-particle 
Liouville equation in Hamilton picture. In the Section 3, in the Liouville picture we substitute the 
canonical distribution function in the Liouville equation for a dissipative system and derive condition 
for non-potential forces. In this section, we consider the Maxwell-Boltzmann distribution function for 
dissipative systems. In Section 4, the example of the canonical distribution for linear oscillator with 
friction is considered. Finally, a short conclusion is given in Section 5. 

2 Canonical distribution from Liouville equation in Hamilton picture 

Let us consider the N-particle classical system in the Hamilton picture. In general, the equation of motion 
the ith particle, where i = 1, N, has the form 

dr^ _ jh dpi _ p 
dt m' dt 
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where is a resulting force, which acts on the iih particle. In the general case, the force is not a potential 
and we can write 

where [7 = U(r) is the potential energy of the system, F- n ^ is the sum of non-potential forces (internal and 
external), which act on the ith particle. For any classical observable A = A(f t ,pt,t), where r = (ri, rjv) 
and p = (pi, ...,pn) in the Hamilton picture, we have 



The Hamiltonian of this system 

N 

H(r.n) = V 

2m 
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« = E£; + c/ « ( 3 ) 



i=l 

is not a constant along the trajectory in the 6N dimensional phase space. From equations (2) and (1) we 
have 



dH _^LpjdU A/ dU p(n)\Pi 
dt ^mdfi df* * /m' 



i.e. 

N 



dt 



i=l 

where i?j = Pi/m. Therefore, the energy change is equal to a power of the non-potential forces P^: 

V(f,p,t) = j2(Ft\v l ). (5) 



i=i 

The N-particle distribution function in the Hamilton picture is normalized using 



J PN (f t ,p t ,t)d N f t d N pt = l. (6) 

The evolution equation of the function pNijtiPtit) is the Liouville equation in the Hamilton picture (for 
the Euler variables) which has the form 

dp N (f t ,p t ,t) _ 

— = -il(r t ,pt,t)p N {r t ,pt,t). (7) 

This equation describes the change of distribution function pjy along the trajectory in 6N-dimensional 
phase space. Here, Q, is defined by 



N d pjn) 



n(*a*) = E->sh (8) 
i=i opi 



and d/dt is a total time derivative (2): 



i q N - a N a 

— = — + yPi— + y Pi—. 

dt dt ~[ m 9r~i ~^ % dpi ' 

If f2 < 0, then the system is called a dissipative system. If CI ^ 0, then the system is a generalized 
dissipative system. In the Liouville picture the function Q is equal to the velocity of the phase volume 
change [10]. 

Let us define a phase density of entropy by 

S(f t ,pt,t) = -k In p N (r t ,p t ,t). 
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This function usually called a Gibbs phase. Equation (7) leads to the equation for the Gibbs phase: 

dS(r t ,p t ,t) 



dt 



kn(f t ,p t ,t). (9) 



Therefore, the function is proportional to the velocity of the phase entropy density (Gibbs phase) 
change. 

Let us assume that the power V(f t ,Pt,t) of the non-potential forces is directly proportional to the 
velocity of the Gibbs phase (phase density of entropy) change ft(r t ,pt,t): 

V(r t ,p t ,t) = kTn(f t ,p t ,t), (10) 

with some coefficient T which is not depend on (ft,pt,t), i.e. dT/dt = 0. 

Using equations (4,5) and (9), assumption (10) can be rewtitten in the form 

dH(r t ,Pt) = T dS(f t ,p t ,t) 
dt dt 

Since coefficient T is constant, we have 

d 



dt 



(H(r t ,p t )-TS(f t ,p t ,t))=0, 



i.e. the value (H — TS) is a constant along the trajectory of the system in 6N-dimensional phase space. 
Let us denote this constant value by T . Then we have 

H(f t ,p t )-TS(f t ,pt,t)=F, 

where dT/dt = 0, i.e. 



In p N (f t ,p t ,t) = H (fu 



Pt 



As the result we have a canonical distribution function 



pN(r t ,p t ,t) = exp^-{T - H{f u pt)) 



kT' 

in the Hamilton picture. The value T is defined by normalization condition (6). 

Note that iV is an arbitrary natural number since we do not use the condition N » 1 or iV — > oo. 

3 Canonical distribution in Liouville picture 

Let us consider the Liouville equation for the N-particle distribution function pjy = pN(r,p,t) in the 
Liouville picture (for the Lagnangian variables): 

„ TV -> o N o 

OPN , Pi OPN , a /- \ 

In general, the forces F{ are non-potential forces. This equation is the equation of continuity for 6N- 
dimensional phase space. Substituting the canonical distribution function 

p N (r, p, t) = exp^ - H(r, p, t 

in equation (11), we get 
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Since Pat is not equal to zero, we have 



^ + yP 1 dH + yp i dH =kT ym 
dt i^i m 1 dpi dpi 



If the Hamiltonian H has the form (3), than this equation leads to 

t— f m v or i / rrf <m 

4 = 1 1 = 1 r 

Substituting equation (1) in this equation, we get the following condition for non-potential forces F-^: 

1=1 1=1 1 

Using notations (5) and (8), we can rewritte this condition in the form 

V(f,p,t) = kTCl(r,p,t). 

As a result we have that the canonical distribution function is a solution of the Liouville equation for 
dissipative and non-Hamiltonian systems if the power of the non-potential forces is proportional to the 
velocity of the phase volume change. 

Let us consider a chain of Bogoliubov equations [11, 6] for the Liouville equation of the dissipative 
systems (11) in approximation 

p 2 {n,pi,r 2 ,p 2 ,t) = p 1 (f 1 ,p 1 ,t)p 1 (r 2 ,p 2 ,t). (12) 

The non-potential forces F^ in equation (1) is a sum of external forces p'f 1 '^ and internal forces p^ n ' l \ 
For example, in the case of binary interactions we have 

N 

Ft ) =F^ e \n,p i ,t)+ £ F^if^r^t). 
In approximation (12) we can define the force 

Fl ~ + Fl + i W 

where 

f[% f (ri , Pi , t) = J df 2 dp 2 pi (f 2 , p 2 , t)F^2 ,l) (n ,pi,f 2 ,p 2 ,t), 

U e ff(ri,pi,t) = J df 2 dp 2 p! (r 2 , p 2 ,t)U(\f 2 - ri|). 
If we consider the 1-particle distribution then Liouville equation (11) in approximation (12) has the form 

dpi Pi dpi d \ 

-m + ^W l + wS Flpl ^ = ^ (13) 

where p\ = p\(f\,p\,t). Let us consider a condition for the non-potential forces 

(pi, F± + FlJ ff ) = mkT — 

In this case, we can derive the 1-particle distribution function (as in Section 3) in the form 

pi(f,p,t) = Aexp-^ f (f- + U{f) + U eff (r)). 
This is a Maxwell-Bolztmann distribution function. 
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4 Canonical distribution for harmonic oscillator with friction 



Let us consider the N-particle system with a linear friction defined by non-potential forces 

i? n) = -IPi, (14) 

where i = 1, N. Note that iV is an arbitrary natural number. Substituting equation (14) into equations 
(5) and (8), we get the power V and the Gibbs phase S7: 



A? 



Condition (10) has the form 

£sl-*t, (15) 

i=i 

i.e. the kinetic energy of the system must be a constant. Note that equation (15) has no friction parameter 
7. Condition (15) is a non-holonomic (non-integrable) constraint [12]. 

Let us consider the N-particle system with friction (14) and non-holonomic constraint (17). The 
equations of motion for this system have the form 

dfi Pi dpi _ dU dG 

W = m' H = -™-W l +X M' (16) 
where the function G is defined by 

N 

2 ' 



1 N 

G(f,P) = 5 (E^ 2 - mkT ) : G (r,P) = 0- (17) 
i=i 



Equations (16) with condition (17) define 6N+1 variables (f,p,X). 

Let us find the Lagrange multiplier A. Substituting equation (17) into equation (16), we get 

dpi , u . flP , 1fi , 
- = -(7-A)^-— . (18) 

Multiplying both sides of equation (18) by pi/m and summing over index i, we obtain 

1=1 1=1 1=1 

Using dG/dt = and substituting equation (15) into equation (19), we get 

= -( 7 -A)W 

r— : J m or,- 

Therefore, the Lagrange multiplier A is equal to 

1 dUs 
J= i j 

As the result, we have the holonomic system which is equivalent to the non-holonomic system (16,17) 
and defined by 

dr\_Pi_ dpi_ _J_vY- ^L\-> 

dt ~ m' dt ~ mkT M Pj '' drJ Pi dr{ 1 Uj 
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Condition (10) or (15) for the classical N-particle system (20) is satisfied. If the time evolution of the 
N-particle system (16) has non-holonomic constraints (17) or the evolution is defined by equation (20), 
then we have the canonical distribution function in the form 



can have a canonical distribution look likes the canonical distribution of the linear harmonic oscillator: 



5 Conclusion 

Dissipative and non-Hamiltonian classical systems can have stationary states that look like canonical 
distribution. The condition for non-potential forces which leads to the canonical distribution function 
for dissipative systems is very simple: thepower of all non-potential forces must be directly proportional 
to the velocity of the Gibbs phase (phase entropy density) change. 

In the papers [13, 14], the quantization of evolution equations for dissipative and non-Hamiltonian 
systems was suggested. Using this quantization it is easy to derive quantum Liouville-von Neumann 
equations for the N-particle matrix density operator of the dissipative quantum system. The condition 
which leads to the canonical matrix density solution of the Liouville-von Neumann equation can be 
generalized for the quantum case by the quantization method suggested in [13, 14]. 
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For example, the N-particle system with the forces 





where 
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